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Abstract. We consider linear and quadratic integrals of motion for general variable quadratic 
Hamiltonians. Fundamental relations between the eigenvalue problem for linear dynamical invariants 
and solutions of the corresponding Cauchy initial value problem for the time-dependent Schrodinger 
equation are emphasized. An eigenfunction expansion of the solution of the initial value problem is 
also found. A nonlinear superposition principle for generalized Ermakov systems is established as a 
result of decomposition of the general quadratic invariant in terms of the linear ones. 



I. An Introduction 

Quantum systems with variable quadratic Hamiltonians are called the generalized harmonic os- 
cillators (see, for example, [2], [13], [T5], [28], [3H], [53]' [SI]) [H5] and references therein). They have 
attracted substantial attention over the years in view of their great importance in many advanced 
quantum problems. Examples include coherent states and uncertainty relations [53], [51], [55], [57] . 
[33], Berry's phase [2], [3J, [28], [39], [62], asymptotic and numerical methods [36], [58], [61], |63j, 
quantization of mechanical systems [22], [23], [23], [25], [33], [35] and Hamiltonian cosmology [3], 
[26] . [27] . [29] . charged particle traps [52] and motion in uniform magnetic fields [6], [10], [15], [37] . 
[33], [33], [3H], [55], molecular spectroscopy [16], [53] and polyatomic molecules in varying external 
fields, crystals through which an electron is passing and exciting the oscillator modes, and other 
interactions of the modes with external fields [25]. Quadratic Hamiltonians have particular appli- 
cations in quantum electrodynamics because the electromagnetic field can be represented as a set 
of forced harmonic oscillators [5] , [15] , [25] , [60] and [73] . Nonlinear oscillators play a central role in 
the novel theory of Bose-Einstein condensation [11] . From a general point of view, the dynamics of 
gases of cooled atoms in a magnetic trap at very low temperatures can be described by an effective 
equation for the condensate wave function known as the Gross-Pitaevskii (or nonlinear Schrodinger) 
equation [3D], [31], [32] and [68]. 

In this Letter we consider the one-dimensional time- dependent Schrodinger equation 

if = Hi P (1.1) 

with general variable quadratic Hamiltonians of the form 

d 

H = a (t) p 2 + b (t) x 2 + c (t) px + d (t) xp, p = -i—, (1.2) 

ox 

where a(t) , b (t) , c (t) , and d (t) are real- valued functions of time, t, only (see, for example, [6], [7J, 
[H], [9], [13], [23J, [23], [25], [38], [50], [59], [76], [77], [78], [79], [83], and [85] for a general approach 
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and some elementary solutions). The corresponding Green function, or Feynman's propagator, can 
be found as follows [5J, [75] : 



V> = G(x,y,t) = 1 e <°oW- 2 +/3 W^+7oW2/ 2 ). 
V2tt^o (t) 



where 

1 ^(t) c{t) 



a (t) 



4a (t) fi (t) 2a (t) ' 



ftW-^l, A(t)=exp( / (r(.s) -,/(*)) r/. 



° w M*W) 2 "(0) A «,(»)) 

and the function /i (t) satisfies the characteristic equation 

At" - r (t) // + 4(7 (t) n = 

with 

a' c /a' c' 

r (i) = h 2c - 2d, a (t) = ab - cd + 

a 2 \ a c 

subject to the initial data 

/i (0) = 0, ^(0) = 2a (0)^0. 

(More details can be found in Refs. [6] and [75].) Then by the superposition principle the solution 
of the Cauchy initial value problem can be presented in an integral form 



1.3 

1.4 
1.5 
1.6 

1.7 



1.9 



^{x,t)= / G(x,y,t) <p(y) dy, Mm ip (x, t) = <p (x) (1.10 

for a suitable initial function (p on R (a rigorous proof is given in [75] and uniqueness is analyzed 
in [5J; another forms of solution are provided by (14.201) and (I7.1ip ). 

A detailed review on dynamical symmetries and quantum integrals of motion for the time- 
dependent Schrodinger equation can be found in [15] and [53] (see also an extensive list of references 
therein). In this Letter, which is a continuation of the recent paper [5J, a natural connection be- 
tween the linear and quadratic integrals of motion for general variable quadratic Hamiltonians is 
established. As a result a nonlinear superposition principle for the corresponding Ermakov systems, 
known as Pinney's solution, is obtained and generalized. We pay also special attention to funda- 
mental relations between the linear dynamical invariants of Dodonov, Malkin, Man'ko and Trifonov 
and solutions of the Cauchy initial value problem (see original works [H], [T5], [53], [SS]). In addi- 
tion this initial value problem is explicitly solved in terms of the quadratic invariant eigenfunction 
expansion, which seems to be missing in the available literature in general. 

The Letter is organized as follows. We start from the standard definitions of the dynamical 
symmetry and quantum integrals of motion, introducing also some elementary tools, in the next 
two sections. Then in Section 4 we describe all linear dynamical invariants for variable quadratic 
Hamiltonians and determine their actions on the solutions of the corresponding time-dependent 
Schrodinger equation. In Section 5 all quadratic quantum integrals of motion are characterized 
in terms of solutions of the generalized Ermakov equation and a detailed proof is given. Their 
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decomposition into products of the linear invariants is derived in Section 6 and explicit actions 
on the solutions are established in Section 7. The last section deals with a related nonlinear 
superposition principle for Ermakov's equations and some computational details are provided in 
Appendices A and B. 



2. Dynamical Symmetry 



In this Letter we elaborate on the following property. 
Lemma 1. // 

4r = i ^' 4r + OH ~ Hf ° = 0> (2 - 1} 

then the function %p l = Oip satisfies the time- dependent Schrddinger equation 

= HH X , (2.2) 

where is the adjoint of the Hamiltonian H. 

When H = H\ this property is taken as a definition of the dynamical symmetry of the time- 
dependent Schrddinger equation ( II. II) (see, for example, [15], [48], [53] and references therein). 
At the same time one has to deal with non-self-adjoint Hamiltonians in the theory of dissipative 
quantum systems (see, for example, [7], [13], [35], [ST], [32] and references therein) or when using 
separation of variables in an accelerating frame of reference for a charged particle moving in an 
uniform variable magnetic field [6J. 

Proof. Partial differentiation 

dip -i d ,. dO , ,„dib , 

'w-'m^-'w^^i (2 - 3) 

= (H ] - OH) 4> + OHip = H^ x 
provides a direct proof. □ 

Definition of the dynamical symmetry is usually given in terms of solutions of the same equation. 
A simple modification helps with the non-self-adjoint quadratic Hamiltonians. 

Lemma 2. The wave functions ip and x, related by 

iP = (e- ti (c ~ d)ds Oj X , (2.4) 

are solutions of the same Schrddinger equation M.1\) - [T7^) , if the operator O satisfies hypothesis of 
Lemma 1. 

Proof. The simplest dynamical invariant, or an operator with the property ( 12. ip . is given by 

O = O (c, d) = e ^ c - d)ds J, (2.5) 

where I = id is the identity operator. (More details are provided in Section 4.) Apply Lemma 1 
twice in the following order 

iP = O (d,c)(O X ) (2.6) 

and use (H^y = H to complete the proof. □ 
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Examples occur throughout the Letter. 



3. Differentiation of Operators and Dynamical Invariants 
Following Lemma 1 we define the time derivative of an operator O as follows 

».% + l t( OH-BV). (3-D 

where W is the adjoint of the Hamiltonian operator H. (This formula is a simple extension of the 
well-known expression p2], [37], [60], [73] to the case of a non-self-adjoint Hamiltonian [7J.) By 
definition 

§-§ + \{O*-*O)- (3,, 

for any dynamical invariant. (In what follows we assume that the invariant O does not involve time 
differentiation.) 

This derivative is a linear operator 

d dOx d0 2 ( q\ 

— (ciOi + c 2 2 ) = ci—— + c 2 — - (3.3) 
at at at 

and the product rule takes the form 

|(0,0 2 ) = 6 ( ^° 2> + j ((dCM H - g' (O.Q,)) (3.4) 

For the general quadratic Hamiltonian, (11.21) . one gets 

^ (0l02 ) = ^ 02 + 0l ^ + (c _ rf)0l02 (3.5) 

and by definition ( 13. 11) 

* ( e a/ <(c- d )d S O Q \ = gcjjMd. ^ 02 + 0l ^ (3 . 6) 

at V / \ at at y 

+ (a + 1) (c - a 1 ) e °^( c - rf ) ds O x 2 . 

If a = —1, we finally obtain 

* (e-JSM* O x 2 ) = e -£(-«0* f ^Io 2 + . (3.7) 

at \ / \ at at J 

This implies that, if the operators 0\ and 2 are dynamical invariants, namely, 

^ + l (O „ tOl) =0, (3,) 

then their modified product, 

E= e -fZ(c-d)d8 Q i0 ^ ( 31 Q) 
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is also a dynamical invariant: 

§ = § + \(EH-^= . (3.11) 

In Section 6 this property allows us to describe connections between linear and quadratic dynamical 
invariants of the time-dependent Hamiltonian (jl.2p . 

4. Linear Integrals of Motion 

All invariants of the form 

dP 

P = A{t)p + B{t)x + C{t), ^ = 0, (4.1) 



(we call them the Dodonov-Malkin-Man'ko-Trifonov invariants; see, for example, [14], [To], [53J, 
and [56J and references therein) for the general variable quadratic Hamiltonian ( U.2p can be found 
in the following fashion. Use of the differentiation formula (13 . 1 j) results in the system [S]: 

A' = 2c(t)A-2a(t)B, (4.2) 

B' = 2b(t)A-2d(t)B, (4.3) 

C = (c(t)-d(t))C. (4.4) 

The last equation is explicitly integrated and elimination of B and A from ( 14. 2 p and ( 14. 3p . respec- 
tively, gives the second-order equations: 

A"- (^ + 2c-2d)A' + 4(ab-cd+~ ^ - £Vj A = °> ( 4 - 5 ) 

B" - (j + 2c-2d)B' + 4 Lb-cd-^(j-^j^B = 0. (4.6) 

The first is simply the characteristic equation ( ll.7p -( |T78l) . (It also coincides with the Ehrenfest 
theorem ([S], [J7]) when c 4-> d.) 

Thus all linear quantum invariants are given by 

P — A(t)p + 2C{t)A 2 { ^ A,(t) x + Co exp Qf (c (s) -d(s)) ds^j , (4.7) 

where A (t) is a general solution of equation (14. 5 j) depending upon two parameters and C is the 
third constant. Study of spectra of the linear dynamical invariants allows one to solve the Cauchy 
initial value problem ([J4], [15], [53] and [56]). 

Theorem 1. (Eigenvalue Problem for the Linear Invariants) If 

/x 2C (t) U (t) - ll' (t) , ln . 

P{t) = u it) p + 2 ^ (t) — x, ( 4 -8) 

i/ien /or any solution A = u(t) of the characteristic equation [T. 7| )-( TQ1) we /iai>e 

P (t) X (x, y,t) = (0) /i (0) A (t) y K (x, y, t) . (4.9) 

The eigenfunctions are bounded solutions of the time- dependent Schrddinger equation U.l\) given 
by 

K (x, y, t) = 1 e <«M- 2 +/3M-/+7(% 2 ) a 10) 
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where X (t) = exp ( J * (c (s) — d (s)) dsj and 

p(t) = 2/i(0)/i («)(a(0)+7 («)), (4.11) 

a(t) = 1 (4.12) 

v 7 4a (t) // (t) 2a (t) v ; 



4 (a (0)+ 7o (*))' 

/?(*) - - (4.13) 
U 2(a(0)+ 7o (*)) 

7 ^ = 7(0)- ,, ^ 2 | 0) my (4.14) 
4 (a (0) +7 (t)) 

(When p, = n with fi (0) = ana 1 (0) = 2a (0) 7^ 0, we obtain the Green function ( fi.^) -( TOp up 
to a simple factor.) 

Proof. The required solution f l4.10p - fl4.14p has been already constructed in [78] and one has to verify 
formula (14.91) only. We have 

PK(x,y,t) = (Ap + Bx)K(x,y,t) (4.15) 
= {{2aA + B)x + (3Ay)K(x,y,t), 

where 

2aA + B = f±£ _ £) A + !£d^ (4.16) 
\2a ji a J 2a 

A (ai _ — 

~ 2a\n A 

provided \i = A. Then j3A — j3 (0) // (0) A and the proof is complete. □ 







Our theorem can be thought as a natural extension to the case of non-self-adjoint variable qua- 
dratic Hamiltonians (11. 2p of a familiar relation between the Green function and linear dynamical 
invariants established in [13], [15], [53] and [56]. The time- dependent factor in the eigenvalue (14. 9p 
corresponds to the statement of Lemma 2. 

In this Letter we are interested, among other things, in a direct verification of Lemma 1 for linear 
and quadratic dynamical invariants. For the general Dodonov-Malkin-Man'ko-Trifonov invariant 
(14 .7p . without any loss of generality, one can separately consider two cases, say, when A(t) = 
with Cq = 1 and when Co = 0. If 

V> a = V eti {c ~ d)d \ (4.17) 

then 

fy± = M e p oic - d)ds + i{c _ d)l p (4 . 18) 
dt dt \ jri v 1 

= (ap 2 + bx 2 + cpx + dxp) ip 1 

+ (c-d)(xp-px)ip 1 = H j ip 1 , (4.19) 
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which takes care of the first case (we have verified once again the statement of Lemma 1 for the 
simplest invariant (I2.5P ). 

In the second case, Cq = 0, we follow Theorem 1 and take a solution A = fi(t) of (14. 5p and 
(11. 7p . which does not have to satisfy the initial conditions of the Green function (II. 9p . Then by the 
superposition principle solution of the corresponding Cauchy initial values problem is given by the 
integral operator 



il>(x,t)= K(x,y,t) X (y) dy, ^(x,0)=/ K {x, y, 0) X (y) dy (4.20) 

J — oo 

with the kernel (|PDjHj£np . Thus 

^ {x, t) = P?P {x, t) = I {PK {x, y, *)) X (y) dy (4.21) 



(we freely interchange differentiation and integration in this Letter; it can be justified for certain 
classes of solutions ([12], jSS], [SB], [S3]))- By choosing (3 (0) \l (0) = 1 in (jUD we obtain 

;>oo 

(x, t) = Pi> (x, t) = e^ c - d ^ ds / K {x, y, t) (y X (y)) dy, (4.22) 



where the second factor obviously satisfies the Schrodinger equation (II. ip (see [78] for details). 
Repeating the first step one completes the proof. Eq. (14.221) provides a spectral decomposion ( [T0J , 
|71j ) of the linear invariant P, which has a continuous spectrum. 

Our consideration shows how Lemma 1 works for the general Dodonov-Malkin-Man'ko-Trifonov 
invariant — application of this invariant to a given solution of the corresponding Cauchy initial 
value problem simply produces another solution with the following initial data: 

V>iOz,0) = P{0)ip{x,0) (4.23) 

/oo 
K(x,y,0) (yx(y)) dy. 
-oo 

The reader can easily connect initial conditions of two solutions (14.201) and (I4.22p with the help of 
an analog of the Fourier transform (see also Ref. [77] for a formal inversion operator in general). 

5. Quadratic Dynamical Invariants 
All quantum quadratic integrals of motion are given by 

dE 

E = A(t)p 2 + B (t) x 2 + C it) (px + xp), — = 0, (5.1) 

(aiL 

and can be found as follows [8]. 



Lemma 3. The quadratic dynamical invariants of the Hamiltonian U.S\) have the form 



E(t) 



c + d) k — k' ^\ 2 | Cq 2 



np-\ x] +—x 

2a ) k 

x exp I / (c — d) ds 



(5.2) 
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where Cq is a constant and function k (t) is a solution of the following nonlinear auxiliary equation: 

\2 



Aab + 



c-d)(c + d) -d -d! 



K = C 



(2a) 



(5.3) 



The structure of quadratic invariants is once again in perfect agreement with Lemma 2. 



Proof. In a general form this result has been established in [8] (see also Refs. [48], [80] and [85] 
for important earlier articles). A somewhat different and more direct proof is presented here. It is 
sufficient to show that the corresponding linear system 

A' + AaC - (3c + d) A 
B' - AbC + (c + 3d) B 
C + 2 (aB - bA) -{c-d)C 
has the following general solution 

(c — d) ds 



0. 
0. 





(5.4) 
(5.5) 
(5.6) 



Ait) 
B(t) 
C(t) 



k 2 exp 



K 



(c + d) K 



2a 



(c + d) k — k' 

k exp 

2a F 



Co 

K 2 



exp I / (c — d) ds 



(c — d) ds 



(5.7) 
(5.8) 
(5.9) 



where Co is a constant and the function k (t) satisfies the nonlinear auxiliary equation (15. 3j) . 

Indeed the substitution 

A(t) = A{t)X{t), Bit) = B(t)X{t), Cit) = C{t)X{t), (5.10) 

where A it) = exp ^ J Q * (c — d) dsj , transforms the original system into a more convenient form 



A' + AaC -2(c + d)A 
B' -AbC + 2(c + d)B 

C' + 2 ( aB - bl) 



0. 
0. 



Letting 



A = k 2 and A' = 2kk' 



in the first equation (15. lip , we obtain 



C 



Then from the third equation (15. 13ft : 



2a 



((c + d) k - K?) 



(5.11) 
(5.12) 

(5.13) 
(5.14) 
(5.15) 



_ e J(c+d)dt ( Ke ~f(c+d)dt\ 



2a dt 
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-k H r 

a 2a at 



2a dt 



f (c+d)dt 



and substitution of fl5.15p - fl5.16p into f)5.12p gives 

d\i 



k 



d_ 

dt 



4abk 2 fi 2 + k^- [ fi [ k 



d 



in the following temporary notations 



Using 
one gets 



v(t) =y(r) 



dt 



- J(c+d)dt 



dt 



+ 8abk 2 fj, ( k 



dfi 
~dt 



k = — e 
2a 



2 / (c+d)dt 



k^ = % q = Aabk 2 
dt dr 



3 4 / (c+d)dt 



d_ 

dr 



qy 2 + 



d ( dy 
dr \ dr 



+ 2,y^ = 

dr 



or 



(see also Ref. 



y (y" + Qy)' + 3y' (y" + qy) = 
for an important special case). By an integrating factor: 

Y T V (y" + w)] = °> y" + iy = ^, 



and a back substitution with the help of 



dr 2 



^3f{c+d)dt 

(2a) 2 



K 



v + 



— c — d I (c + d) — d — d' J k 



results in the required first integral of the system: 



d 
dl 



(2a) 



K 



k' + ( 4a6 + ( c-d) (c + d) - d - d' ) k 



which gives our auxiliary equation (15. 3p . 

The last term in (I5.16P can be transformed as follows 



1 d 
2adt 



J(c+d)dt 



k d 



-2 f(c+d)dt 



-2 J (c+d)dt 



2a dt 

d 



- f(c+d)dt 



dy 



dr V^dr 



dy 



Tr ] ~ qy + ~ 



-2f(c+d)dt (yy" + (y> 



Co 



K 



(c + d) K, 



2a 



Co 



K 



in view of (I5.18p -f l5.19p and ( 15. 22ft . We have also utilized a convenient identity 



dy 



dr 2a \ dt 



dK 



— (c + d) k ) e 



J(c+d)dt 
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Thus 

V- (c + d)K\ 2 C 
2a J K- 

and the proof is complete. □ 



The case a = 1/2, b = u> 2 (t) /2 and c = d = corresponds to the familiar Ermakov-Lewis- 
Riesenfeld invariant [2TJ, [13], [H], [13], [IB]; see also Refs. [IB] and |17] . (The corresponding 
classical invariant in general is discussed in Refs. [SD] and [85].) Many examples of completely 
integrable generalized harmonic oscillators and their integrals of motion are given in Ref. [3J. 

For a positive constant, Co > 0, the quantum dynamical invariant (15. 2p can be presented in the 
standard harmonic oscillator form [7], [8], [70], [71] : 



2 

where 

u; (t) = exp 



^(o^at^ + at^aCt)), (5.27) 



(J^(c-d)ds), w = 2 v / C^>0, (5.28) 

3 «> = i^-^^r")-^ < 5 -»> 



/ ,/Un K 1 — (c+d)K,\ K 



2k 2ay/uJo J y/UJodx' 

and k is a real- valued solution of the nonlinear auxiliary equation (I5.3p . Here the time-dependent 
annihilation a (t) and creation a) (t) operators satisfy the canonical commutation relation: 

a{t)a ] (t) - a ] (t)a(t) = 1. (5.31) 

The oscillator-type spectrum of the dynamical invariant E can be obtained now in a standard way 
by using the Heisenberg-Weyl algebra of the rasing and lowering operators (a "second quantization" 
], the Fock states): 

a (t) * n (x, t) = Vn * n _i (x, t) , a f (t) * n (x, t) = Vn + 1 ^n+i (a, *) , (5.32) 
£ (t) * n (x, t) = u (t) ( n + - J ^ n (x, t) . (5.33) 



2, 

The corresponding orthonormal time-dependent eigenfunctions are given by 

k' — (c + d) k 2 



* n (x, t) = exp ( i ^— ) w n (£) (5.34) 



in terms of Hermite polynomials 

^ = C n e-?/ 2 H n (0 , £ = ex, (5.35) 

where 



|C„| 2 = -^-r, e 2 = ^ = (5.36) 
1 1 ^2 n nV 2k 2 k 2 v ; 

Their relation with the original Cauchy initial value problem is discussed in Section 7 (see The- 
orem 2). In addition the n-dimensional oscillator wave functions form a basis of the irreducible 
unitary representation of the Lie algebra of the noncompact group SU (1,1) corresponding to the 
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discrete positive series V J + (see [59], [HI] and [75]). Operators f l5.29f> — ( T5.30D allow us to extend these 
group-theoretical properties to the general quadratic dynamical invariant (15.27p . 

When Co = 0, the dynamical invariant (I5.2p is given as the square of a linear invariant up to a 
simple factor, which resembles the case of a free particle. If Co < 0, one deals with the Hamiltonian 
of a linear repulsive "oscillator" . 



6. Relation Between Linear and Quadratic Invariants 



By Lemma 3 the operators p 2 , x 2 and px + xp form a basis for all quadratic invariants. Here we 
take two linearly independent solutions, say fi 1 = A\ and fi 2 = A 2 , of equations ( 11. 7ft and ( 14. 5 ft and 
consider the corresponding Dodonov-Malkin-Man'ko-Trifonov invariants ( 14. 7ft : 

P 1 = A 1 p + B 1 x, P 2 = A 2 p + B 2 x. (6.1) 

Introducing the following quadratic invariants 

E x = Pi e -ti {c - d)d \ E 2 = P 2 e -fo{c-d)d^ E 3 = (PxP 2 + P 2 P X ) e ~ ^ c ~ d)ds (6.2) 

with the help of (I3.10P as another basis, one gets 

E = C\E\ + C 2 E 2 + C 3 E 3 (6.3) 

= [C x Pl + C 2 P 2 2 + C 3 {P X P 2 + P 2 Px)] exp (-J q (c- d ) ds ^) 
for some constants C\, C 2 and C3. As a result the following operator identity holds 



(c + d) k — k' 

k p H x 

y 2a 



. C 2 



exp 
,2 



(c — d) ds 



[d (A lP + Six) 2 + C 2 (A 2 p + B 2 xf 
+C 3 {{Axp + Bxx) {A 2 p + B 2 x) + (A 2 p + Sax) {A ± p + B x x))\ 



where 



Thus we obtain 



Ax = f-ix, Bx 



x ox]) ( — / (c — d) (is 
2c/i x — n'x 



2a 



An 



M 2 ! 



^9 



2c/x 2 - ^ 2 
2a 



k 2 = (Cifil + C 2 /i 2 + 2C 3 /i 1 /i 2 ) exp ^—2 y (c — d) ds 



(6.4) 



(6.5) 



(6.6) 



as a relation between solutions of the nonlinear auxiliary equation (I5.3P and the linear characteristic 
equation ( |4.5p . In addition the substitution 



fix = Kx exp (c — d) ds*j , /x 2 = ^ 2 exp ^y (c — d) ds^j 



(6.7) 



transforms the characteristic equation (14. 5 p into our auxiliary equation ( 15. 3 p with Co = 0. Finally 
a general solution of the nonlinear equation is given by the following "operator law of cosines" : 

k 2 (t) = Cxk\ (t) + C 2 k\ (t) + 2C 3 Kx (t) K 2 (t) (6.8) 
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in terms of two linearly independent solutions K\ and k 2 of the homogeneous equation. The constant 
Cq is related to the Wronskian of two linearly independent solutions K\ and k 2 : 

v 2 



dC 2 - Cl = C -J^ -, W («!, k 2 ) = /d/c'a - Ki« 2 (6.9) 

(more details are given in Appendix A). This is a well-known nonlinear superposition property of 
the so-called Ermakov systems (see, for example, [10], [H], [21], [10], [H], @2], [15], [51], [15], [67], 
[69] . [71] and references therein). Here we have obtained this "nonlinear superposition principle" 
(or Pinney's solution) in an operator form by multiplication and addition of the linear dynamical 
invariants together with an independent characterization of all quantum quadratic invariants, which 
seems to be missing, in general, in the available literature (see also |43j and [15] for an important 
classical case). An extension is given in the last section. 

It is worth noting, in conclusion, that the linear invariants of Dodonov, Malkin, Man'ko and 
Trifonov ([H], [15], [S3] and [5B]) can be presented as follows 



P 1 = U lP + ( - C + d \ K a 1 ^ xjexplj (c-d)ds), (6.10) 

P 2 = [k 2P + (C + ^ ~ exp QT (c - d) ds^j (6.11) 

in terms of two linearly independent solutions, K\ and k 2 , of the homogeneous equation ( 15. 31) when 
Co = 0. Comparing these expressions with the form of the quadratic invariant ( 15.21) at Cq = 
(no "potential", a "free particle"), when the operator square is complete, one can treat the linear 
invariants as "operator square roots" [12] of the special quadratic invariants (see also Lemma 2 
regarding a convenient common factor). 



7. Quadratic Invariants and Cauchy Initial Value Problem 

Our decomposition (16.31) of the quantum quadratic invariant in terms of products of the linear ones 
not only results in the Pinney solution (I6.8l) - (l6.9p of the corresponding generalized Ermakov system 
(15. 3p in a form of an "operator law of cosines" , but also provides a somewhat better understanding, 
with the help of Lemma 1 and properties of the linear invariants discussed in Section 4, how the 
quadratic invariants act on solutions of the original time-dependent Schrodinger equation. Indeed 
by ( 14.201) for two different solutions, say A\ = fi l and A 2 = /i 2 , of the characteristic equation, ( 14. 5ft . 
we have in an operator form: 

^ (x, t) = K x (t) [ Xl (y)\ = K 2 (t) [ X2 (y)\ (7.1) 

in view of uniqueness of the Cauchy initial value problem (see, for example, Refs. [8] and [78]). 
Then by (T47221 : 

= e Io(^)ds Ki {yXi) ) = e ffc-d)ds K2 ( yx j (7 2) 

and 

= e -^ c - d)ds [C 1 P^ + C 2 P 2 2 ^ + C 3 (P 1 P 2 + P 2 P 1 )^P] (7.3) 
= eJ?(~0* [ ClKl (y 2 Xl ) + C 2 K 2 {y\ 2 )\ 
+C 3 [P 1 (K 2 (y X2 ))+P2(K 1 (y Xl ))}. 



<p 1 = Kr 1 (0) [K 2 (0) (y Xa )] 
p 2 = K 2 l (0) [K x (0) (y Xl )} 



DYNAMICAL INVARIANTS 

We have 

K 2 (VX2) = K x M , 
K\ (yxi) = K 2 (<p 2 ) , 

with the help of an analog of the Fourier transform. Therefore 

Pi (K 2 (y X2 )) = e'° (c - d)ds K x ( y(Pl ) , P 2 (K x (y Xl )) = K 2 (yip 2 

and finally 

Eij) = e^~ d ^ [CxKi (y 2 Xl ) + C 2 K 2 (y 2 X2 ) + C, (K, ( m ) + K 2 (y<p 2 ))] , 
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J 2J > 



(7.4) 



where each term satisfies the corresponding Schrodinger equation. By the superposition principle 
we arrive at a new solution in a complete agreement with our Lemma 1. The corresponding initial 
data follow from (17.41) at t = and the time-evolution operator (11.101) can be applied (see also 
Eq. (I7.17P for an eigenfunction expansion). 



On the second hand one can always expand a square-integrable solution (I4.20p of the Cauchy 
initial value problem in the standard form 



ip{x,t) 



K(x,y,t) X (y) dy 

) 

^c n {t) V n (x,t), 



(7.5) 



n=0 



where 



(x, t) ip (x, t) dx 



(7.6) 



(we use the asterisk for complex conjugate) by the Riesz-Fisher theorem [70], [71], [72] due to 
completeness of the eigenfunctions (I5.34I) - (I5.36I) at all times. Then by the Fubini theorem: 



x{y) 



I ®* n (x,t) K (x,y,t) dx) dy 

J — oo / 



(7.7) 



and the second integral can be evaluated with the help of Eqs. (14.101) . (I5.34p (I5.36j) and (IB. 1H as 
follows 



/oo 
%(x,t)K(x,y,t) dx 
-oo 

/ 



(7.8) 



jn e -i(n+l/2)<p 

( ( 



\ 



i/i 



x exp 



V 



7 + 



n I K.l\ 2 I ( ti\K'—K',K\ 

\°o {-) + y — Ta — ) J 



exp ( - / (d — c) ds 



o 

\ \ 

2 



V 



n i ki \ 2 i ( k\k'-k' x k\ 

°0 {-) + y Ta ) J 



y 



J 
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/ 



x exp 



/3 2 (0) «? (0) VC~ t 



n ( Kl\ 2 I ( KiK'-k[k\ 2 

V C ° + (, 2a ) J 



H, 



P(0) K\ (0)C 1/4 i/ 



\ 



Here «i = /zexp ( f Q * (d — c) dsj and k are the corresponding solutions of auxiliary equation (15. 3p 
with C = and C ^ 0, respectively, with the Wronskian W ( 

^1) ^0 — K\K K-^K and 

K K\K K-^K 



tan ip 



+ 



«l 2ay / Co 

2 



K~^K K^ 

2a 



constant 



(7.9) 



(7.10) 



by ()A.4j) . (The computational details are left to the reader; our equation (IT. lOj) is equivalent to the 
classical Ermakov invariant [21].) We may choose /3 (0) K\ (0) = 1 and arrive at the following result. 

Theorem 2. (Eigenf unction Expansions) Solution of the Cauchy initial value problem (i4.20\) 
in L 2 (R) can be obtained as an infinite series of multiples of the quadratic invariant eigenf unctions 



V> (x, t) = J2°n (*) *n *) 



n=0 



where the time- dependent coefficients are given by 

5 ^ 1/2 



c n (t) 



e -i(n+l/2)<p exp _ 



V / vr2 n n! 

POO 

x / exp(z^ 2 )e- 5 V/2 Fn( ^ )x(y) dy 



(d — c) (is 



(7.11) 



(7.12) 



Here K\{t)and n(t) are real-valued solutions of the homogeneous and nonhomogeneous auxiliary 
equations A5.3\) , respectively, with the Wronskian W (t) — K\K K-^K. The phases (p (t) and 7 (t) are 
determined in terms of these solutions as follows 



(p = arctan 



K K^K ^\ ^ 

K\ 2aJ~C~o 



dt k 2 ' 



^ + - = 
dt k\ 



and 



a 



1/4 



6 2 



> 0, 



^ = 7 + 



K \ K KlKjK 



(7.13) 
(7.14) 

(7.15) 
(7.16) 



2^/Co 2aK\ 

are constants. A spectral decomposition of the quadratic invariant E in the space of L 2 -solutions is 
given by 



E (t) t/j (x, t)=u (t) ^ (*) (n + |^ * n (ar, t) 



(7.17) 
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where the "frequency" u (t) is defined by A5. 28\) . 
It is worth noting that 



d_ 

It 



2^/C^K 1 \2aJ_ 



(7.18) 



by (I7.14p and ( 1A.6I) . which means that the phase factor (17.161) in front of y 2 in the second integral 
of Eq. (17.121) is indeed a constant. The second equation (17.1 3p follows from (17.1 8p with the help of 
( ESP and dZH. 

Finally by choosing in Eqs. fl7.lip - fl7.12p a special (square) integrable initial data of the form 

Xm (y) = exp (-z£y 2 ) e^ y2l2 H m (8y) (7.19) 

we conclude in view of the orthogonality property of Hermite polynomials that the time-dependent 
wave functions 

^ m (x,t) = D m exp Qjf (d-c) ds^j (7.20) 
Xe -i(m+i/2Mt) ^ m ( X)t ) 
are particular solutions of the Schrodinger equation fll.ll) - fll.2p for arbitrary constants D m : 

i*bjl = H1> mi Ei> m = u(m+'^ji> m . (7.21) 

They are also eigenfunctions of the quadratic invariant E corresponding to the discrete "spectrum" : 

(E)m = J_ r m E^ m dx=\D m \ 2 u Jo (m + ^j, (7.22) 

see (I5.33p . (It can be verified by direct substitution; the details are left to the reader.) The explicit 
wave functions f )7.20p are derived here without separation of the variables with the aid of our 
Theorem 2 and certain variants of the Ermakov invariant. If \D rn \ = 1, solution (17. lip takes the 
form 

°° /"OO 

1> (x, t) = Yl i n e-^ n+1 ' 2 ^ ^ [x, t) / r n (V, 0) X (V) dy (7.23) 

in terms of the time-dependent wave functions ( 17.20P provided that 

(c(0) +<i (0)) K (0)-.-(0) 
V ; k(0) 4o(0)«(0) V 1 

(see also [85] for the path- integral method). The traditional operator approach (for the parametric 
oscillator) is presented in [48] and/or elsewhere (see also [I] and [3"9"]). 



8. A General Nonlinear Superposition Principle for Ermakov's Equations 



The Pinney superposition formula fl6.8p - fl6.9p allows one to construct solutions of the nonlinear 
auxiliary equation (15.31) in terms of given solutions of the corresponding linear equation. In general 
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we take two different solutions, say K\ and k 2 , of the generalized Ermakov equations (15.31) with 



a 



(1 2) 

0/0 and consider two quadratic invariants: 

c + d) Ki 



Ei(t) 



Ki p 



2a 



x 



+ 



a 



(i) ' 

x 2 



E 2 (t) = 
Their arbitrary linear combination 



(c + d) K 2 — K-'o 
K 2 P~\ X 



2a 



+ 



a 



(2) ' 

x 2 



X(t) 
X(t) 



(8.2) 



(8.3) 



D l E 1 (t) + D 2 E 2 (t) = E (t) 

(D\ and D 2 are constants), is also a quadratic invariant given by (15.21) for a certain solution k of 
the nonlinear auxiliary equation (15.31) . Thus the following operator identity holds 

2 r,„ _ 

U) 



(c + d) k — k' 
k p H X 



D 1 



+D 2 



2a 



KiP + 



, 2 
+ x 

K 



[c + d) Hi 



2a 



x 



+ 



a 



(i) ■ 

x 2 



(c + d) K 2 — k' 2 

K 2 V H 7s x 

2a 



+ 



a 



(2) 

x 2 



K 2 



and in a similar fashion we arrive at a general nonlinear superposition principle for the solutions of 
Ermakov 's equations (15. 31) : 

k 2 (t) = D x k\ (t) + D 2 k 2 2 (t) , (8.5) 

where the constant Co is given by 



Co - C^Dl - C^D\ = D 1 D 2 



( 2 ) n2 



(2a) 



r (l) K 2 , r (2) K l 
°0 . 2 ' i 



(8.6) 



with k 2 ) being the Wronskian of two solutions (see Appendix A). One can also derive this 

property by adding the corresponding solutions (I5.7p -( [5T91) of the original linear system (I5.4j) ( f5T6T) 
or with the help of the Pinney formula (I6.8p - (l6.9p . The details are left to the reader. 
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Dodonov, Vladimir I. Man'ko, Svetlana Roudenko and Kurt Bernardo Wolf for support, valuable 
comments and encouragement. I am grateful to Peter G. L. Leach for careful reading of the manu- 
script — his numerous suggestions have helped to improve the presentation. 



Appendix A. Required Transformations 

Equating coefficients in front of the operators p 2 , x 2 and px + xp in (I6.4p - fl6.5p with the help of 
(16. 7p one gets 

k 2 = C x k\ + C 2 k\ + 2C 3 kik 2 , (A.l) 



(c + d) K — k' 
2ci 



n 2 



Co 

K 2 



(A.2) 
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(c + d) K\ — K± 

2a 



C 2 



(c + d) K 2 — «2 

2a 



+2C7; 



[(c + d) ACi — k'J [(c + d) k 2 — k' 2 ] 



(c + d) K — k' 
2a 

„ (c + d) Ki — k[ _ (c + d) k 2 — k' 2 
— Gi«i h C 2 K 2 — 



(A.3) 



+C 3 



2a 

(c + d) k 2 - k' 2 
2a 



+ k 2 - 



2a 

c + d) K\ — /Cj 
2a 



respectively. Multiply (IA.2[) by (lA.ip and use flA.3[) in the left hand side in order to obtain 
a result of elementary but rather tedious calculations. 



as 



A general nonlinear superposition principle for Ermakov's equations, given by (I8.5p - (l8.6p . is 
derived from the quadratic invariant identity (j8.4p in a similar fashion. In addition the constant in 
the right hand side of (18. 6p (see also (17.101) ) can be verified by direct differentiation as follows 



d 
dt 



Ta) +C ° UJ + ° W 

W'--w)+2W -d®± 
a J \ K i K 



(A.4) 



(2a)" 
2W 
(2a7 



W--W 

a 



2W 
J2af 



t i i \i i i i \ 

[K\K 2 ~ K 1 K 2J (KiK 2 — K-iK 2 ) 

~ a 







with the help of auxiliary equations (I5.3p . (It is a natural generalization of the classical Ermakov 
invariant [21].) En route, we derive an identity: 



wY 



,(2) Kl 



2a \2a J 

which may be considered as an extension of the familiar Abel theorem, occurring when C { 

(2) 

Cq = 0, to the nonlinear Ermakov equations. Another identity, 



(A.5) 



(i) 



d 


"«2 / 


'WY 


dt 


Ki \ 


>A 




2a 


' (W 




'-^ 


\2a 



(A.6) 



.(1) / K 2 



«1 



^1 -C7-I-] + cf * 



«2 



has been useful in Section 7 with C, 



(i> 



0, c 



(2) 



Cq and k 2 = k when the right hand side reduces 



to a multiple of the Ermakov invariant (see Eq. (I7.18P ). 
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Appendix B. An Integral Evaluation 
In Section 7 (see Eq. (17.81) ) we use the following integral 

/oo 
e -#(*-v) a Hn ( ay) dy 
-oo 

X n+l{ ) \ (A 2 - a 2 ) J 

which is equivalent to Eq. (30) on page 195 of Vol. 2 of Ref. [19] 
polynomials), or Eq. (17) on page 290 of Vol. 2 of Ref. [20] . 
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